Medium polarization and finite size effects on the superfluidity of the inner crust of 

neutron stars 
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The J 5o pairing gap A associated with the inner crust of a neutron star is calculated, taking into 
account the coexistence of the nuclear lattice with the sea of free neutrons (finite size effects) , as well 
as medium polarization effects associated with the exchange of density and spin fluctuations. Both 
effects are found to be important and to lead to an overall quenching of the pairing gap. This result, 
whose quantitative value is dependent on the effective interaction used to generate the single-particle 
levels, is a consequence of the balance between the attractive (repulsive) induced interaction arising 
from the exchange of density (spin) modes, balance which in turn is influenced by the presence of 
the protons and depends on the single-particle structure of the system. 
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I. INTRODUCTION 

Neutron stars are possible remnants of supernova, an 
explosion signaling the death of a massive star when it 
has run out of nuclear fuel, displaying a spectacular rapid 
increase in brightness. 

Theoretical models consistent with the experimental 
findings testify to the fact that a neutron star has a thin 
atmosphere and three internal regions referred to as : 
the outer crust, the inner crust, the core. In the outer 
crust, matter consists of spherical atomic nuclei and elec- 
trons. At the bottom of this region, neutrons begin to 
drip out of nuclei, thereby producing a neutron gas be- 
tween nuclei. It is generally accepted that in the den- 
sity range O.OOlpo < P < 0.5po, where p = 2.8 • 10 14 
g cm corresponds to nuclear saturation density, neu- 
trons stars display a superfluid inner crust, consisting of 
a Coulomb lattice of nuclei permeated by a gas of free 
neutrons. The superfluidity of the inner crust has im- 
portant consequences concerning different aspects of the 
physics of neutron stars, such as its heat capacity and 
thus its cooling rate. It is also responsible for a number 
of macroscopic quantum phenomena, such as quantized 
vortices (in rotating neutron stars) and quantized mag- 
netic flux tubes (in magnetized neutron star cores). The 
interaction of vortices with nuclei forming the Coulomb 
lattice in the inner crust is thought to be connected with 
the presence of glitches in the pulsation periods of neu- 
tron stars and the postglitch relaxation of the associated 
period (cf. e.g. ref. [l[). 

For a consistent description of these phenomena it is 
necessary to take into account finite size effects arising 
from the coexistence of finite nuclei and free neutrons. 
Moreover, polarization processes can renormalize the in- 
teraction between particles in an important way, and thus 
strongly influence the superfluid properties of the system. 



For example, in uniform neutron matter a strong quench- 
ing of the So pairing gap arising from the polarization 
of the medium is generally predicted 0], although no 
consensus has been reached concerning the actual inten- 
sity of this effect. Hence, superfluid properties can be, 
in principle, correctly described only if medium polar- 
ization processes (e.g. self-energy, induced interaction, 
vertex corrections) are properly taken into account @ . 

The present paper represents a step along the path 
which leads to this ambitious goal. We divide the inner 
crust into spherical Wigner-Seitz (WS) cells of different 
densities, and take into account both finite size effects 
(arising from the presence of a nucleus at the center of 
the cell) and, among all possible medium polarization 
effects, those associated with the interaction induced by 
the exchange of medium fluctuations, which are expected 
to give the largest contribution at the low densities typ- 
ical of the system under consideration. We shall use a 
theoretical scheme which can be applied both to finite 
nuclei [1] and to the inner crust of a neutron star, the 
only essential difference between these two systems be- 
ing the value of the Fermi energy. A preliminary study 
along these lines, limited to a sjngle value of the neutron 
density, was published in ref. ||. 



II. THE MEAN FIELD 

We consider five regions of different densities in the in- 
ner crust previously studied by Negele and Vautherin in 
the WS cell approximation f(|. The number of protons, 
Z, and the radius of the cell, RwSi have been taken from 
their study, and are listed in Table |TJ We determine the 
single-particle states inside each cell in the Hartree-Fock 
(HF) approximation, making use of either the SLy4 or 
the SkM* Skyrme interactions @, H . These HF calcu- 
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lations are analogous to those previously performed by 
other groups [t| [l(| ■ We use wavefunctions vanishing at 
the edge of the cell, working on a grid with a mesh size 
of 0.1 fm. The boundary condition affects the single- 
particle density only in a region of w 3-4 fm from the 
edge of the cell (see Fig. [1]). Single-particle levels up to 
100 MeV have been computed. We have also performed 
calculations without protons, choosing the Fermi energy 
Ep so as to obtain about the same number of neutrons 
in the cell. We have verified that in the case where no 
protons are considered we reproduce, to a good accuracy, 
the pairing gap found in infinite uniform neutron matter 
PH . The resulting densities and the neutron numbers ob- 
tained in the various cells without protons are reported in 
Table [Til We also report the value of the effective mass, 
m e ff, associated with the effective interaction used. The 
neutron and proton density distributions calculated with 
the SLy4 interactions turn out to be quite similar to those 
calculated in ref. [6| . They are shown in Fig. Q] in the 
case p n = 0.01 fm~ 3 . 



SLy4 



SkM* 



z 


N 


Rws 
(fm) 


40 


458 


42.2 


50 


1048 


35.8 


50 


1314 


33.2 


50 


1740 


27.6 


40 


1474 


19.6 


40 


458 


42.2 


50 


1314 


35.8 


50 


1738 


33.2 



TABLE I: Proton number Z, neutron number N and radius 
Rws for the WS cells containing a finite nucleus. 
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FIG. 1: Neutron and proton mean field densities calculated 
with the SLy4 interaction in the cells 1364 Sn and 1314 X (where 
X stands for "no protons") . 



TABLE II: Neutron density p n (and its ratio with the neutron 
saturation density po), Fermi momentum Uf, Fermi energy 
Ef, effective mass m e f f /mt ar e, neutron number TV and radius 
Rws for the WS cells without protons. 



III. COLLECTIVE MODES 

We have computed the excitation spectrum of the sys- 
tem within the framework of the Random Phase Approx- 
imation (RPA) in configuration space using the particle- 
hole basis states Ijpj^ 1 , JM). The residual interaction 
was determined as the second derivative of the mean 
field energy functional, that is, it was derived from the 
same two-body Skyrme force which determines the mean 
field, except for the two-body Coulomb and the two-body 
spin-orbit interactions which have not been included in 
the calculations. RPA phonons are classified according 
to their total angular momentum J and parity 7r, while 
an index k = 1,2,... runs over the states having the 
same J n values with increasing energy. Because of the 
spin-orbit interaction, each phonon \k, J n ) is actually an 
admixture of S = and S = 1 states. However, since the 
spin-orbit effect is very small for free neutrons, which 
dominate the RPA response, phonons with non-natural 
parity are almost pure S = 1 states, while natural parity 
states are either almost pure S = or almost pure 5=1 
states. 

To obtain convergence in the matrix elements of the 
induced interaction, one needs to consider angular mo- 
menta up to about J = 30. For a given value of J, the 
matrix elements converge rather rapidly with the phonon 
energy, which enters directly into the denominators of 
such matrix elements, (cf. Eqs. §5§ below). The maxi- 
mum size of the particle-hole space is 4000 particle-hole 



3 



configurations. We have checked in few cases that this is 
enough to achieve convergence in the calculation of the 
pairing gap. 

The HF+RPA calculation without protons could be 
tested against a calculation performed in neutron mat- 
ter. We have checked that in our configuration space we 
are able to reproduce the response of the infinite uni- 
form neutron matter to the external fields V ex t(r) — e lq ' r 
(5 = channel) and V ext (r) — e lq ' r a z (5=1 channel). 
Analytic formulae for such responses have been given in 
ref. [l2j • The response function in the WS cells has been 
evaluated according to 



S(q,E) 



n=0 



(n|e i «- f e s |0)| 2 L(S,£; n 



(1) 



where \n) are the (mean field or RPA) excited states of 
the system, L is a Lorentzian function used to smooth 
out the computed discrete response, and S = 1 or a z 
for 5 = and 5=1 respectively. Because of the spher- 
ical symmetry adopted, we have expanded the external 
operator in multipoles and obtained the total response 
function adding the contributions from each multipolar- 
ity For this comparison the SLy4 parametrization has 
been used. The analytical response function for uniform 
neutron matter is compared with the WS cell approxi- 
mation in Figs. ([2^) and (J^b) for a neutron gas with a 
Fermi momentum fc^=0.39 fm _1 and A;f=1.16 fm _1 re- 
spectively. These two systems correspond to the lowest 
and highest density regions of the inner crust of a neutron 
star (WS cells with 508 and 1480 neutrons respectively). 
The comparison has been performed for a transferred mo- 
mentum q equal to 0.51 fm _1 for the lower density region 
and 0.76 fm _1 for the higher density region. 

The overall agreement found between the numerical 
and analytic results testifies to the fact that the WS cell 
approximation is a reliable tool to calculate the response 
of the free neutron gas which permeates the nuclear lat- 
tice. We assume that the same is also true for the case 
of a WS cell with a nucleus in the center. Within this 
context, a more detailed study of the RPA response will 
be given in a different paper. 



IV. INDUCED INTERACTION MATRIX 
ELEMENTS 

The matrix elements of the induced pairing interaction 
in the 1 5q channel 
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FIG. 2: The response functions S per unit volume as a func- 
tion of the energy, computed in the WS cells at &f = 0.39 
fm _1 (N = 508) (a) and k F = 1.16 frn -1 (N = 1480) (b), are 
compared to the analytical results in uniform neutron matter. 
Calculations are made with SLy4 interaction. Both unper- 
turbed (HF) and RPA responses are compared. Solid lines 
correspond to the uniform neutron system. The linear mo- 
mentum transferred from the external field is q — 0.51 fm~ 
(a) and q = 0.76 fm" 1 (b). 



where |a6;00) (a = n a l a j a ), indicates a normalized two- 
particle state coupled to total angular momentum zero, 
are computed in Bloch-Horowitz perturbation theory [l3|, 
EI 01 El: 

(ab;00\v lnd \cd;00) = 

(ab; 00\H coup i\int) (int\H coup i\cd; 00) 



Eq — Ei n t 



(3) 



where Ei nt is the energy of the intermediate state made 
by two-quasiparticles and one phonon, 



Es,, 



Eb + Ed + E\. 



(4) 



The quasiparticle energies are denoted by Eb and Ed, 
while E\ is the energy of the exchanged phonon and Eq is 
the binding energy of a Cooper pair. In keeping with the 
results of previous studies (cf. e.g. 0,[I3j[3[l3) we have 
set Eq equal to — 2Ap, where Ap represents the value 
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of the state-dependent pairing gap obtained solving the 
HFB equation, averaged over the single-particle states 
lying in the energy range Ep±2 MeV. We also computed 
the second time ordering (see Eq. @), associated with 
the intermediate state of energy Ei nt = E a + E c + E\. 

Hcoupi represents the particle-vibration coupling 
Hamiltonian. It reads 



H, 



coupl 



mi A 



fmi Mmj^A + A mi O } 



(5) 



with 



' mi 




E 

nj 



^mjin^-nj 



u mnij 1 nj 



). 



(6) 

where A^ mi creates the particle- hole excitation \m(i)~ l ), 
0\ is the phonon creator operator, X^j and are the 
forward and backward amplitudes of the phonon |A) asso- 
ciated with the particle- hole component v m jin 
is the RPA antisymmetrized matrix element of the resid- 
ual interaction v res . We treat the latter in the Landau 
approximation. Thus 



= [F (r)+F^r)f-r'}5(r-^) 
+ [G (r) + G' (r)T-f'}a-a'S(r- 



(7) 



Exchange terms of v res have not been considered. More- 
over, as far as the computation of the induced interaction 
matrix elements is concerned, the terms associated with 
the momentum dependent terms of the NN effective in- 
teraction have been neglected [ij. An analysis of the 
influence of these terms will be presented in future work. 
We shall only consider the t z -t' z term, in keeping with the 
fact that we are interested in the neutron-neutron pairing 
interaction. In uniform neutron matter, RPA phonons 
are characterized by their spin: density (S = 0) modes 
are associated uniquely with the spin-independent part 
of the interaction, given by Cf{t) = F (r) + F^r), while 
spin (5 = 1) modes are produced by the spin-dependent 
part, C G (r) = G (r) + G' (r). In the WS cell, one has 
to take into account also the proton degree of freedom 



-1) contribution to the RPA response. More- 



over, because of the spin-orbit interaction, each phonon 
is an admixture of S = and S = 1 states. While cou- 
pling to non-natural parity modes can take place only 
through the spin-dependent part of the interaction, cou- 
pling to natural parity modes can receive contributions 
from both channels, although the spin- independent part 
is usually the dominant one. 

The explicit expressions for the vertices / and g in the 
S = and S = 1 channels, coupling the neutron single- 
particle states to the K-phonon of angular momentum 
and parity J 7r , are 




K J"M , S=0 



= (fm'\H coup i\jm, kJ v M) 

= (j'm'\i J YjM\jm) J dr^(r)[(F + F^)Sp% n (r) + 
(Fo-F^)Sp^ p (r)]Mr) 

(8) 




KfM , S=1 



(j'm'\H coup i\jm, KJ n M) 

.7+1 

]T (j'm'\i L \Y L ((l) x a] JM \jm) 

L=J-1 

dr^,(r)[(Go + G' )Sp^ Ln (r) + 
(G ~G' )5p^ Lp (r)]Mr), 

(9) 

(f) v {r)/r is the radial single-particle wavefunction associ- 
ated with the quantum numbers v = {nlj}. In the above 
expressions Spj^ n and 6pj^ p are the neutron and proton 
contributions to the transition densities associated with 
the RPA modes \nJ v ). For the 5 = channel one finds 



ph 



Mr)/r Vp\\i J Yj\\j h ) 



r x 



V2J+ 1 



(10) 



while in the 5 = 1 channel we have 



ph ^ ' 

W)/r (j P \\i L [Y L (n)xa]j\\j h ) 



V2J+ 1 



(11) 

The neutron and proton contributions are obtained sum- 
ming over the neutron or proton particle-hole states. 
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Since the sign of the denominator in the expression of the 
induced pairing matrix element is always negative (see 
eq. [3|), and since for scattering vertices there is an ad- 
ditional negative sign in the S — 1 channel as compared 
to the S = channel [4j, it follows that the exchange of 
vibrations gives rise to an attractive pairing interaction 
in the S = channel, and to a repulsive interaction in 
the S = 1 channel. Hence, the induced interaction aris- 
ing from the exchange of normal parity modes involves 
in general a mixture of attractive and repulsive contri- 
butions, while for non- normal parity modes only repul- 
sive terms are present. One can show, however, that the 
repulsive contribution to the diagonal matrix elements 
(nlj, nlj; 00\vi n d\nlj, nlj; 00) associated with normal par- 
ity modes vanishes. These diagonal matrix elements are 
shown in Fig. [S^a) for the case where protons are con- 
sidered ( 1364 Sn), and in Fig. g^b) for the case without 
protons, ( 1314 X), corresponding to the density p n = 0.01 
fm~ 3 . They have been calculated with the SLy4 force. 
The matrix elements of the two systems are of the same 
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order of magnitude («0.01 MeV), except for a few, much 
larger matrix elements which appear when protons are 
present. The large matrix elements connect two-particle 
states based on resonant single-particle states which have 
positive energy e n y and are spatially localized within 
the nucleus (particularly on the nuclear surface). The 
largest contributions to these matrix elements arise from 
collective phonons with a well defined surface charac- 
ter, both for density and spin modes. As an example, 
let us focus on the diagonal matrix elements associated 
with the single-particle state (n — 3, 1 = 7,j = 15/2; 
e 3,7.i5/2 = 2.02 MeV), whose values are equal to -0.49 
MeV for normal parity, and to +0.31 MeV for non-normal 
parity (see arrows in Fig. [3Ja)). The transition densities 
of the phonons which give rise to the most important 
contribution to each of these large matrix elements are 
shown in Fig. [4l together with the single-particle wave- 
function 03,7,15/2- 
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FIG. 4: (a) Proton and neutron transition densities asso- 
ciated with the = 2 + normal parity phonon with energy 
4.46 MeV for the cell 1364 Sn. The single-particle wavefunction 
4>nij = ^3,7,15/2 is a l so shown, (b) The same, for the J n = 15 + 
(L=14) non-normal parity phonon at 15.27 MeV. There is 
no proton contribution to this transition density within our 
particle-hole space. 



FIG. 3: Induced pairing interaction diagonal matrix elements 
(in the BCS limit) for the cells 1364 Sn (a) and 1314 X (b), com- 
puted with the SLy4 effective interaction. Positive (negative) 
matrix elements correspond to abnormal parity (normal par- 
ity) modes. Note the different scale of the two figures. The 
arrows indicate the two matrix elements discussed in the text, 
and associated with the transition densities shown in Fig. H 



The matrix elements obtained making use of the SkM* 
interaction are qualitatively similar, but there are signif- 
icant quantitative differences. A schematic plot of the 
diagonal induced interaction matrix elements obtained 
with the two Skyrme forces used in the present paper 
is shown in Fig. [5] It can be seen that the matrix ele- 
ments associated with the SkM* force are larger, in ab- 
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solute value, than those associated with the Sly4 force. 
Moreover, the spin modes have matrix elements larger 
(in absolute value) than the density modes, while they 
are about the same for SLy4. This can be understood 
considering the Landau parameters factors Cf and Cq 
as a function of density in uniform neutron matter. As 
it can be seen in Fig. O in particular at low density, the 
parameters associated with the SkM* have much larger 
values (in absolute value); moreover Cq is larger than 
Cp, while the opposite is true for SLy4. 




FIG. 5: Schematic representation of the BCS diagonal in- 
duced interaction matrix elements for the two cells 1314 X (left 
panels) and 1364 Sn (right panels) obtained making use of SLy4 
(upper panels) and SkM* (lower panels) effective interactions. 
The grey area represents the range of the computed matrix 
elements (panel (a) corresponds to Fig. ([3b), while panel (b) 
corresponds to Fig. (|3^,) except for the few, larger and scat- 
tered matrix elements). 



V. THE PAIRING GAP 

The pairing matrix elements of the bare and induced 
interaction were used in a HFB calculation to study 
the superfluid properties of the system in the 1 So chan- 
nel. The HFB calculation is self-consistent only in the 
particle-particle channel, while single-particle eigenstates 
and eigenvalues are kept fixed and equal to those com- 
puted in the HF calculation used to produce the mean 
field discussed previously. 

We first consider pairing correlations at the mean field 
level, employing the Gogny interaction Q as the bare 
interaction Vbare- This is rather well justified as long as 
the neutron density is smaller than, or of the order of 
0.01 fm~ 3 (kp = 0.7 fm _1 ), because the pairing gaps ob- 
tained with the Gogny or the bare interaction are quite 
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FIG. 6: Diagonal matrix elements of the Gogny interaction 
at the Fermi level (solid line) compared to the residual inter- 
action form factors Cf (dashed lines) and Co (dash-dotted 
lines) for the uniform neutron matter when the SLy4 (a) and 
the SkM* (b) Skyrme interactions are employed, (cf. Eq. 0) 



similar in that density range. Instead for higher neutron 
densities and in the nuclear interior the Gogny interac- 
tion produces somewhat larger gaps. While this differ- 
ence is important for a quantitative study of the pairing 
gap both in infinite and in finite systems [21], it may be 
neglected in the present context, where a qualitative dis- 
cussion of the role pairing interaction plays in the neutron 
star crust is the main issue. 

The state dependent pairing gap A nn nj associated 
with the two-particle state \nlj, n'lj) is given by 



A 



"5 ^.V^l? 17 " 

n 2 nLl 2 n 2 < \ 1 



yi 

2/2 32 n' 2 l 2 j 2 



n 2 n 2 l 2 j 2 ' V 1 

(n 1 l 1 j 1 ,n'l 1 j 1 ;00\v tot \n 2 l2j2,n' 2 l 2 j2; 00), 

(12) 

where f7*y and V£ t ■ are the quasiparticle amplitudes ob- 
tained solving the HFB equations in the pairing channel, 
and v to t = Vbare + v ind- We have taken into account the 
single-particle levels contained in an energy window go- 
ing from Ep ± 9 MeV for the lowest density WS cell to 
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Ep ± 20 MeV for the highest density WS cell. Increasing 
the window by 1 MeV or including one more multipolar- 
ity changes the pairing gap by less than 0.01 MeV in all 
cases. 



A. Cells without protons 

The values of A;? obtained with and without the in- 
duced interaction on top of the Gogny interaction, and 
calculated in the WS cells without protons, are shown 
by filled and open triangles in Fig. [71 It is seen that in 
uniform neutron matter the induced interaction quenches 
the value of the pairing gap, and that the suppression is 
stronger at large densities. 
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FIG. 7: Pairing gap calculated at the Fermi energy for the WS 
cells without protons employing the SLy4 (a) and the SkM* 
(b) Skyrme interactions. Empty and filled triangles refer to 
the results of the microscopic calculation with the Gogny in- 
teraction only and with the Gogny plus the induced interac- 
tions respectively. The solid line corresponds to a BCS calcu- 
lation of the pairing gap in neutron matter, while the dashed 
line corresponds to a simpler calculation (see Appendix A) 
which employs a Landau approximation of the residual inter- 
action both for RPA and for induced matrix elements. 

When the SLy4 interaction is used, the induced inter- 



action reduces the maximum value of the pairing field 
by 25%, bringing it from about 3.2 MeV to about 2.4 
MeV. Moreover the gap goes rapidly to zero around kp w 
1.2 fm , compared to the progressive decrease occurring 
when only the Gogny interaction is acting. The effects of 
the induced interaction associated with the SkM* case, 
are even more drastic. In this case the pairing gap goes 
to zero already at kp ss 0.9 fm -1 and the gap shows a 
maximum value of about 1.2 MeV. 

The results reported above can be checked and bet- 
ter understood, calculating the pairing gap in infinite 
uniform neutron matter with and without the induced 
interaction using the simple expressions collected in Ap- 
pendix A, which lead to the curves shown in Fig. \7\ The 
agreement with the much more involved, discrete calcu- 
lation performed in the WS cells is quite satisfactory, 
also considering that some approximations have been in- 
troduced in the expressions used in infinite matter. The 
stronger screening effect obtained with the SkM* interac- 
tion is related to the different balance between spin and 
density modes already discussed in the previous section 
in relation with Fig. [S] 



B. Cells with protons 

We now turn our attention to the calculation including 
protons. The pairing field obtained for all the computed 
cells is shown in Fig. [8l The circles refer to the calcula- 
tion with the nucleus at the center of the cell, while the 
triangles refer to the same gaps without nucleus already 
shown in Fig. [JJ 

In Fig. Ufa) we show the results obtained with the 
SLy4 interaction. One can see that the presence of the 
nucleus reduces the pairing gap by about 100-200 keV, 
both with and without the induced interaction. The re- 
sult without the induced interaction is similar to that 
obtained in ref. [HI, where it was shown that the pair- 
ing gap is reduced by about 200 keV close to the Fermi 
energy. This effect was attributed to the fact that pair- 
ing is quenched inside the nuclear volume, where the local 
Fermi momentum is larger than in the outer neutron gas. 
This can be seen in a clear way in Fig. [9fa) where we 
show the spatial dependence of the pairing gap for the 
case with nucleus at p n =0.01 fm -3 . This is obtained 
from our quantum calculation as explained in ref. [22| : 
first we construct the non local pairing field A(ri,r2); 
then we perform a Fourier transform repect to the rela- 
tive coordinate r\ — r*2, obtaining the gap as a function 
of the relative momentum k and of the center of mass co- 
ordinate Rcm, A(fc,i?cAf) (we neglect the weak depen- 
dence on the angle between k and Rcm)', finally we define 
a local pairing gap Ai oc (R C m) = A(k F (RcM), Rcm), 
where kp(RcM) is the local Fermi momentum. From 
this figure one observes that in fact the local pairing gap 
is reduced from 2.7 MeV far away from the nucleus, to 2.4 
MeV in the interior, and thus the mean value of the gap 
in the whole cell is also reduced (the reduction would be 
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stronger using the bare V14 Argonne interaction instead 
of the Gogny interaction, cf. Figs. 5 and 10 of ref. [ill]). 




FIG. 8: Pairing gap calculated at the Fermi energy for the 
different WS cells. Empty and filled triangles refer to the 
calculation without protons, performed with the Gogny in- 
teraction only (empty triangles) and with the Gogny plus the 
induced interaction (filled triangles), already shown in Fig. 
Empty and filled circles correspond to calculation with 
protons (nucleus placed at the center of the cell), performed 
with the Gogny interaction only (empty circles) and with the 
Gogny plus the induced interaction (filled circles) . Mean field 
and RPA modes are computed making use of the SLy4 (a) 
and of the SkM* (b) interactions. 

Turning now to the calculations including the induced 
interaction, considering the results obtained in neutron 
matter, where the pairing gap goes to zero for kp > 1.2 
fin -1 , one might expect that the pairing would almost 
vanish inside the nucleus (cf. crosses in Fig. (Hb)). In- 
stead the actual local pairing gap resulting from the mi- 
croscopic calculation (solid curve in Fig. [9jb)) shows 
that the gap inside the nucleus has a value of about 1.3 
MeV, which represents a much smaller reduction. This 
is to be interpreted as a result of the enhancement of the 
attractive (negative) induced matrix elements due to the 
exchange of nuclear collective surface modes, as it was 
discussed in the previous section (cf. also Fig. [3]). The 
gap in the outer neutron gas remains larger than in the 




10 15 20 

R rM (fm) 



FIG. 9: Pairing field A(k F ,R CM ) for the cell 1364 Sn, calcu- 
lated with the SLy4 mean field and the Gogny pairing inter- 
action (a), with the SLy4 mean field and the Gogny+induced 
pairing interaction (b), with the SkM* mean field and the 
Gogny pairing interaction (c) and with the SkM* mean field 
and the Gogny+induced pairing interaction (d). The crosses 
in (b) and (d) are obtained dividing the pairing field shown 
in (a) and (c) respectively, by the quenching factor calculated 
in uniform neutron matter. 



interior of the nucleus, so that the resulting average gap 
in the cell is somewhat smaller than in the case without 
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the nucleus. 

Considering now the results obtained with the SkM* 
interaction, one can observe two main differences com- 
paring Figs. H](a) and (b): first, as already discussed, the 
screening effect due to the interaction is much stronger 
than in the SLy4 case; secondly, the effect of the nucleus 
now leads to an increase of Ap when the induced interac- 
tion is included. This is due to the fact that, as we have 
seen above, the SkM* interaction leads to a much larger 
screening compared to SLy4 in uniform matter at low 
density. In fact, the induced interaction leads to rather 
similar values of the gap inside the nucleus for the two 
forces, while spin modes reduce the gap in the external 
neutron gas much more with the SkM* than with the 
SLy4 force (cf. Fig. O^b) and (d)). In summary, the 
action of the induced interaction acts in a similar way 
for SkM* and SLy4 in the interior of the nucleus, where 
the action of spin modes (which at saturation density is 
about the same for the two forces, cf. Fig. is strongly 
counteracted by surface modes; while in the outer, low- 
density neutron gas the induced interaction associated 
with SkM* suppresses the gap much more effectively than 
in the SLy4 case. 



VI. CONCLUSIONS 

The effect of the induced interaction arising from the 
exchange of density and spin modes of the system results 
in a quenching of the neutron pairing gap, both in the 
case of uniform neutron matter and in the case of the 
inner crust of a neutron star. 

The pairing gap obtained is found to be strongly de- 
pendent on the effective nucleon-nucleon interaction em- 
ployed to determine the single-particle properties of the 
system and its linear response (we performed the calcula- 
tion making use of SLy4 and SkM* Skyrme interactions). 
A much stronger suppression of the gap is obtained in 
uniform neutron matter when the SkM* interaction is 
used, due to the stronger residual interaction associated 
with this NN effective force and to the different balance 
between spin and density modes. 

When a nuclear cluster is present in the cell, the local 
Fermi momentum increases and the pairing gap calcu- 
lated with the bare interaction is reduced in the interior 
of the nucleus. The induced interaction quenches the gap 
both inside in the nucleus and in the outer neutron sea: 
the relative amount depends on the interaction, but the 
reduction inside the nucleus is much smaller, than that 
occurring in neutron matter at the corresponding Fermi 
momentum. 
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VIII. APPENDIX A 

The detailed calculations of the pairing gap performed 
in the WS cell without protons can be checked against a 
much simpler estimate in neutron matter, using a similar 
assumption for the residual interaction, which is taken 
into account within the Landau approximation, neglect- 
ing the momentum dependence of the Landau parame- 
ters. However, the response function will also be esti- 
mated within the Landau approximation, while in the 
microscopic calculation in the WS cell the full RPA re- 
sponse has been considered. The matrix elements of the 
induced interaction in infinite neutron matter are ob- 
tained summing the contributions vfj~P and vf^ 1 in the 
5 = and S = 1 channels. They are obtained multi- 
plying the square of the Landau parameters C F (density 
modes) and Cq (spin modes) by the integral of the RPA 
response function R s=0 (q) (density modes) and R s=1 (q) 
(spin modes) [23|, [24[ ■ For density modes one has 



C 2 F {k F ) 
k\k 2 



dqqR s =o(q), (13) 



where k\ and k 2 are the momenta of the states exchang- 
ing the vibrations. The integration limits are equal to 
q-min = \k\ — k 2 \ and to q max — k\ + k 2 . One has 



1 N(0)L(g) 

2 1 - C F (k F )L{q) 



(14) 



where L{q) is the Lindhardt function and N(0) is the 
density of single-particle states at the Fermi energy, 
N{0) = %™2 ■ The contribution from spin modes is cor- 
respondingly obtained from R s=1 , where 



1 N(0)L(g) 
21 -C G (k F )L(qY 



(15) 



and from 



3C%(k F ) r> 
k\k 2 



dqqR s=1 (q), (16) 



where the factor 3 is associated with the three spin pro- 
jections in the S = 1 channel. 
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The pairing gap is then calculated solving the BCS 
equations with the complete interaction, v to t(ki, k%) — 
Vbare + vf r Td + v fnd • ^ n Fig- El we have compared these 



simple expressions with the microscopic calculations in 
the Wigner cell without nuclei. 
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